Quantum interference in attosecond transient absorption 
of laser-dressed helium atoms 
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We calculate the transient absorption of an isolated attosecond pulse by helium atoms subject to a 
delayed infrared (ir) laser pulse. With the central frequency of the broad attosecond spectrum near 
the ionization threshold, the absorption spectrum is strongly modulated at the sub-IR-cycle level. 
Given that the absorption spectrum results from a time-integrated measurement, we investigate the 
extent to which the delay-dependence of the absorption yields information about the attosecond 
dynamics of the atom-field energy exchange. We find two configurations in which this is possible. 
The first involves multi photon transitions between bound states that result in interference between 
different excitation pathways. The other involves the modification of the bound state absorption 
lines by the IR field, which we find can result in a sub-cycle time dependence only when ionization 
limits the duration of the strong field interaction. 

PACS numbers: 32.80.Qk,32.80.Rm 
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Attosecond transient absorption (ata) studies provide 
a way to push our understanding of the energy transfer 
between electromagnetic fields and matter to the sub- 
femtosecond time scale p], Q . They are an all-optical at- 
tosecond metrology that complements methods based on 
the measurement of charged particles, such as attosec- 
ond streaking 0, Q and electron interferometry Q. 
Like those methods, high time-resolution is gained in an 
ATA spectrometer by using attosecond extreme ultravi- 
olet (xuv) pulses that are synchronized to the field os- 
cillations of an infrared (ir) laser pulse. The first ATA 
experiments used an attosecond pulse to probe a valence 
electron wave packet (ewp) created by ionizing an atom 
with a strong IR laser pulse 0, [?J . Recent ATA experi- 
ments have used attosecond pulses as a pump that cre- 
ates an ewp which is then probed by a moderately strong 
IR field. In this xuv-pump/lR-probe configuration the 
frequency-resolved transient absorption signal varies as a 
function of the IR intensity, duration, and the sub-cycle 
timing between the two fields [8l4Io|. These experiments 
raise the possibility of studying time-dependent absorp- 
tion down to the attosecond time scale. 

In this theoretical study, similar to recent experiments 
[ToL [Tl| , we consider helium atoms that are excited by an 
isolated attosecond pulse (iap) with a central frequency 
near the ionization threshold, together with a delayed 
few-cycle ir pulse. We elucidate the key features of the 
resulting delay-dependent attosecond absorption (ddaa) 
spectra. These features derive from the fact that the 
IAP "starts the clock" by exciting the system at a well- 
defined time, and that it makes an ewp that is essentially 
independent of the pump-probe delay. The DDAA mea- 
surement is a spectrogram that records the interference 
between different excitation pathways which lead to the 
same absorption or emission processes. Because the at- 
tosecond pulse is locked to the ir field oscillations, it is 
suggestive that many of the features in the spectrum are 
modulated at half the IR laser period (Tir), about 1.3 fs. 
However, the absorption spectrum results from the time- 



integrated response of everything that happens after the 
initial xuv excitation, which complicates the analysis. 
Our main concern in this paper is the connection between 
these oscillations and the real-time attosecond dynamics. 

We discuss two distinct manifestations of attosecond 
dynamics in the DDAA spectrum. First, near the ioniza- 
tion threshold the delay-dependent absorption exhibits 
fully modulated interference fringes at half the IR pe- 
riod. We show these are due to the interference between 
two pathways, separated in time, which give rise to the 
same dipole resp onse. This "which way" quantum path 
interference [12j can be used to time-resolve two pho- 
ton transitions between excited states, and shows how 
the absorption at xuv frequencies is altered on a sub- 
cycle time scale. We find Tir/ 2 oscillations with a sim- 
ilar origin in absorption features associated with non- 
linear XUV+IR processes, which appear as light-induced 
structures in the ATA spectrum. These fringes illustrate 
how energy can be exchanged between excitation modes 
(dressed states) that exist only when the XUV and IR 
fields overlap in time. The second feature results from 
the IR driven sub-cycle AC Stark shift of the lower-lying 
2p and 3p resonances [TO]. The time-dependent Stark 
shift leads to a dispersive line-shape, which results from 
interference between the time-dependent dipole induced 
before, during, and after the IR pulse. We show that, 
due to the integrated nature of the absorption spectrum, 
a measurement of the instantaneous energy shift is only 
possible when ionization limits the interaction time. 

The essential features of DDAA spectra can be under- 
stood at the single atom level. The energy lost or gained 
by a light field in the interaction with an atom can 
be described by a frequency-dependent response func- 
tion S(uj,td) = 21m[d(uj)E*(uj)], where td is the pump- 
probe delay (we use atomic units unless otherwise indi- 
cated). d(uj) and £(oS) are the Fourier transforms of the 
time-dependent dipole moment d(t) and the total driv- 
ing field E(t). The dipole moment is obtained by solving 
the time-dependent Schrodinger equation (tdse) in the 
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FIG. 1: Single atom response function S(uj,td) in helium 
where td is the time delay in IR cycles between the IR laser 
pulse (800 nm, 3 x 10 12 W/cm 2 , 4 cycles, cos 2 envelope, sine- 
like carrier envelope) and the attosecond pulse (330 as, cen- 
tered at 25 eV). The IR intensity oscillations are shown in 
black in the top panel. 



single-active-electron approximation [13j . 

Fig. 1 shows a typical response function vs pump-probe 
delay. The 25 eV, 330 as IAP pump pulse has a band- 
width of 5.5 eV, which means that it overlaps all the 
singly excited and low-energy continuum states of the 
He atom. The FWHM of the 800 nm probe pulse is 11 
fs, corresponding to 4 optical cycles (Tir = 2.7 fs). The 
IR field is of moderate intensity, so that it can not it- 
self excite the atom in its ground state. The pump and 
probe fields have parallel linear polarizations. We include 
a dipole dephasing time T 2 = 65 fs in the calculations by 
smoothly windowing the dipole moment. We have ver- 
ified that using a longer T 2 does not change any of our 
conclusions. Positive (negative) values of S(u,td) mean 
that a dilute gas will absorb (emit) energy at frequency 
uj as the dipole-driven source term in the wave equation 
will be out of (in) phase with S(uj) fl3l. [lij. 

We begin by briefly describing the main features in 
Fig. 1. For large positive delays, td > 4 Tir, the IAP 
arrives after the end of IR pulse and the absorption spec- 
trum exhibits only one-XUV-photon transitions from the 
ground state to the np states. When the two pulses 
overlap, —4 Tir < td < 4 Tir, these absorption features 
are strongly modified. In this same delay range, light- 
induced structures (LISs) appear, for example, between 
the 2p and 3p resonant lines. These are associated with 
two-photon (XUV±IR) processes that transfer population 
from the ground state to non-dipole coupled ("dark") s 
and d states, and have recently been observed experi- 
mentally [11]. Finally, when the IAP arrives before the 
IR pulse, td < —4 Tir, the dipole established by the xuv 
pulse undergoes free decay until the IR pulse arrives and 
strongly modified the dipole, which induces sidebands on 
the main resonance features. 

Many of the features in S(uj,td) show a modulation 
at one half the laser period. Of particular interest are 



fringes that are present both when the pulses overlap 
and when the IAP arrives before the IR. An example is 
seen in Fig. 1 between 24.2 and 24.8 eV near the label 
u 2p + 2cj" , also shown in more detail in Fig. Efb). These 
fringes are caused by quantum interference between two 
distinct pathways for establishing the same coherence be- 
tween the ground state and a group of up states (n > 5), 
near threshold. The process is diagrammed in Fig. [2fa). 
The direct pathway is an XUV-driven lst-order process 
that populates all of the np states simultaneously in a 
manner that is independent of the delay. The indirect 
pathway is a 3rd-order process in which amplitude in the 
2p states is transferred to an np state by a two-lR-photon 
process. The role of the indirect pathway via the 2p state 
is confirmed by a test calculation which dynamically elim- 
inates the 2p state during the time-propagation of the 
tdse, in which these fringes disappear. The two inter- 
fering processes, since they are driven by the xuv field 
and the IR field respectively, happen at different times 
depending on the value of td- The resulting interference 
fringes are analogous to those observed in photoelectron 
spectra [6[, however, in DDAA the multiphoton process 
need not result in ionization. 

Having identified the interfering pathways it is 
straightforward to write down the constructive interfer- 
ence condition as a function of delay: 



(cj np - uJ2 P )\t d - t \ + Acj) = 2tt/c, 



(1) 



where uo np and U02p are the energies of the np and 2p 
states respectively, and k is a positive integer ;6]. Two 
additional parameters enter: Acj) is a phase due to the 
two-lR-photon transition and to is the time when the two 
photon transition probability peaks. It is reasonable to 
choose to = — Tir/4 due to the sine-like IR carrier wave 
used in the calculation. This formula correctly predicts 
a decreasing slope for the fring es as \td — to\ increases, 
corresponding to increasing k in Eq. (pQ). The fit to the 
interference fringes using Acj) = 0, one set of integers k 
that start with k = 1 at to ? an d no other free param- 
eters, is excellent for td < 0, as shown in Fig. [21(b). 
For td > the biggest contribution to the indirect path- 
way comes from whatever IR field maximum follows the 
xuv pulse (never more than one half cycle away). The 
slope of the fringes should therefore saturate for td > 0, 
which is indeed what we observe in Fig. 1. We can thus 
conclude that the strength of the two photon transition 
peaks at the local intensity maxima of the IR field. This 
supports a time-dependent picture of multi-photon ab- 
sorption below threshold as a process that follows the 
sub-cycle oscillations of the IR electric field. 

One application of this 1st vs. 3rd order interference 
process is illustrated in Fig. [2fa). Realizing that the in- 
direct pathway is enhanced due to a near one-lR-photon 
resonance with the intermediate 3s state, it seems pos- 
sible that we can use this to determine Acj) in Eq. (pQ). 
For an 800 nm pulse we can safely assume that Acj) = 
because the intermediate 3s state lies above the energy 
reached by absorbing one IR photon from the 2p state. 
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FIG. 2: (a) Diagram of direct (one-photon) and indirect 
(three-photon) pathways to the same final mp (or ri2p) states. 
Comparison of absorption maxima and the constructive inter- 
ference condition (white dashed curves) for a narrow range of 
frequencies and IR wavelengths of (b) 800 nm (same as Fig. 1) 
and (c) 660 nm. 



In Fig. [2fc) we have repeated the calculation in Fig. 1 
using a wavelength of 660 nm. Since the absorption of 
one photon now reaches an energy above the 3s state, we 
expect this phase should shift by tt. In Fig. [2fc) we use 
= 7r for the same to and achieve an excellent fit to 
the calculation. The interference fringes can therefore be 
used to measure the phase of an intermediate resonance if 
the pulse wavelength can be scanned [l5|. Another, more 
speculative application, is as a probe of decoherence. A 
short decoherence time would reduce the fringe contrast 
for large negative delays, from which the time scale for 
decoherence could be inferred directly. 

We can apply a similar analysis to the LISs as well. 
In the overlap region (— 4Tir < td < 4Tir) each LIS is 
associated with an s or d bound state which is coupled to 
the ground state by one XUV photon and the absorption 
or emission of an additional IR photon. We use a plus or 
minus label for LIS that are one photon above or below 
the final s or d state. LISs associated with coupling of the 
ground state to the 3s, 3d, and 2s states can be seen near 
21.3 eV (3s-), 21.7 eV (3d_), and 22.1 eV (2s+), respec- 
tively. Similar LISs are visible above threshold (3s+, 3d+ 
and 4s+). The half cycle modulations in the LISs come 
from the 2co>ir coupling of plus and minus LISs, which 
have an atomic state as an intermediate resonance. For 
example, the 3d- is two-photon resonantly coupled to 
the 3d+ (visible near 25 eV around td = —ITjr), yield- 
ing oscillations that are are exactly out of phase with 
each other. These couplings can also be confirmed by 
examining the tilt of the LIS interference fringes. Eq. (1) 
predicts that the fringes of the 3d- have a negative slope 
with respect to delay, which is indeed the case in Fig. 
1. However, the fringes in 2s + LIS have a positive slope, 
characteristic of being above the final 2s state. Fits to 
Eq. (1) using the LIS energies show that the two photon 
transitions between LISs also follow the IR electric field 



oscillation. 

A time-dependent perspective on these oscillations is 
provided by the observation that the coupling between 
the one-xuv-photon resonant rip state and nearby dark 
states is in a regime where the rotating wave approxima- 
tion breaks down, and that the short duration of the IAP 
yields a well-defined phase of the IR field at excitation. 
In such situations it has been demonstrated that the final 
state populations are sensitive to the phase of the laser 
at the time of excitation [l6|. Here we have shown that 
the DDAA spectrum shows the same sensitivity. 

Another source of fast dynamics in the transient ab- 
sorption spectrum is the sub-cycle AC Stark shift of the 
bound state energies [10]. Figs. 3(a) and (b) show the 
evolution of the 2p and 3p line shapes with delay, which 
include both emission and absorption. Qualitatively, this 
results from the redistribution of XUV energy across the 
resonances due to the perturbation of the dipole by the IR 
field. Similar dispersive line shapes have been discussed 
previously in connection with the control of exciton po- 
larizations of a quantum dot using picosecond laser pulses 

In a time-domain picture, the dispersive shape of the 
rip absorption lines is caused by multiple contributions 
to the dipole moment: (i) a perturbed response during 
the IR pulse, and (ii) free decay of the coherence after 
the IR pulse ends [H, [20|. The relative importance of 
these terms is determined by the IR field intensity and 
delay, as well as T2. If the state amplitude is unchanged 
by the IR field then the main contribution to the dipole 
moment comes from the free decay which, however, ac- 
quires a phase shift, #s, equal to the integrated optical 
Stark shift during the time when the IR pulse acts. In 
this approximation the response is given by 

S(u, t d ) « £(uj, T 2 ) [cos(0s) + (A • T 2 ) sin(0 s )] (2) 

where A is the difference uo — uj np and C(uo,T2) is the 
Lorentzian line shape in the absence of the IR field. For 
positive delays 1 < td < 4Tir we find that the line 
shape is well described by this simple form. We see 
that the shift of the absorption is in the direction of the 
Stark shift, but the magnitude of the shift is propor- 
tional to 6s /T2 and is therefore constrained by the natu- 
ral linewidth. Thus, for these delays, the 2p and 3p lines 
have Stark shifts of opposite sign, as shown in Fig. 3(c), 
but the magnitude of these shifts does not generally equal 
the instantaneous Stark shift. 

At delays near zero (— 4Tir < td < 1 ^ir) the line- 
shape becomes much wider than I/T2 and more complex 
than the simple dispersive shape discussed above. Here 
depopulation of the excited states by the IR field plays 
a larger role (see Fig. 3(d-e)), and the response is dom- 
inated by the perturbed dipole decay. For the 3p state 
the loss is predominantly by one photon ionization to 
the continuum and the dipole response exists for only a 
few laser cycles after it is excited. This rapid ionization- 
induced dephasing allows us to probe the real-time sub- 
cycle dynamics using a delay-dependent method: during 
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FIG. 3: Response as a function of frequency at around 2p (a) 
and 3p (b) resonance (lineouts of Fig. 1 at different delays). 
Energy shifts retrieved from the maximum response (c) for 
2p (black) and 3p (red) states. Half-cycle modulation of ab- 
sorption strength shown in integration of response function 
(black and red), as well as in final populations at the end of 
the pulses (green and blue) for 2p (d) and 3p (e) sates. 

this time (near -2Tir in Fig. 3(c)) we observe a large 
transient shift of the 3p feature of about 0.19 eV, compa- 
rable to the ponderomotive energy of 0.18 eV, and strong 
half cycle oscillations in the peak position, showing the 
rapid variation in the energy of the resonance as it fol- 
lows the IR field oscillations. These observations are in 
line with other investigations which show that the Stark 
shift is difficult to observe optically in the absence of 
rapid dephasing [7[. This means that observing fast dy- 
namics via ionization-induced dephasing requires careful 
control of the IR intensity. In the continuum, however, 
all states have a very short dephasing time, so we always 
observe an overall ponderomotive shift of the absorption, 
though the response is very weak. 

The behavior of the peak of the 2p lineshape is more 
complex than the 3p, in part because the 2p is never 
fully depopulated during the IR pulse (Fig. 3 (d)), and in 
part because of the strong coherent coupling between the 
2p state and the nearby dark states. The 2p lineshape 
therefore always results from both perturbed and free 
decay. The sub-cycle oscillations of the 2p position are 
most likely due to the which-way interference with the 
np states discussed in the context of Fig. 2, judging from 
the negative and changing slope of the fringes between 
— 3Tir and 3Tir (as seen in Fig. 1). As such they do not 
contain information about the sub-cycle AC Stark shift. 

For delays td < — 4Tir the IAP and IR act separately: 



the dipole established by the IAP oscillates freely until 
the IR pulse strongly perturbs it, greatly altering the 
dipole amplitude and phase in a few IR cycles. This per- 
turbed free polarization decay has been observed previ- 
ously at optical frequencies pjj [21|. It yields sidebands, 
seen clearly above the 3p and below the 2p lines, with a 
characteristic hyperbolic shape that depends on the delay 
as 1/td and shows no attosecond time scale dynamics. 

To conclude, we have identified two configurations 
where attosecond dynamics can be observed in attosec- 
ond transient absorption spectra: which-way interference 
and the sub-cycle AC Stark shift. Since experimental 
observations always result from propagation in a macro- 
scopic gas, we also solve the Maxwell wave equation 
(mwe) for the time propagation of light fields through 
the atomic helium gas medium, in which the polarization 
and ionization source terms are obtained by solving the 
single-atom tdse [13]. The absorption/emission proba- 
bility after propagation of a dilute gas is almost identical 
to the single atom result in Fig. 1, exhibiting both the 
half-cycle modulations, the LISs, and the emission fea- 
tures discussed above. This indicates that our DDAA pre- 
dictions can be observed using current attosecond tech- 
nology, especially as it regards half-lR-cycle oscillations. 

We have demonstrated that interference fringes in 
DDAA spectra result from the coherent addition of two 
quantum paths that lead to the same dipole excitation. 
They reveal the time delay between the initial excita- 
tion and a later, iR-field driven multiphoton transition. 
This was found to be true both for transitions between 
bound states as well as between excitation modes (the 
LISs) that are observable only when the pulses overlap. 
We expect therefore that they will be a general feature 
in more complex systems [22| , and could be observed be- 
tween resonant states embedded in a continuum, as long 
as the lifetimes are longer that the IR period. Using this 
argument in reverse, the interference fringes visible in 
the DDAA spectrum could be used as a precise timing de- 
vice, a probe of decoherence, or a phase meter when an 
intermediate resonance is involved. 
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